a=[5,3,1/3;1,1,1];
b=[100;100];
disp(pinv(a)*b);